Piezoresistance in silicon at uniaxial compressive stresses up to 3 GPa 
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The room-temperature longitudinal piezoresistance of n-type and p-type crystalline silicon along 
selected crystal axes is investigated under uniaxial compressive stresses up to 3 GPa. While the 
conductance (G) of n-type silicon eventually saturates at « 45% of its zero- stress value (Go) in 
accordance with the charge transfer model, in p-type material G/Go increases above a predicted 
limit of « 4.5 without any significant saturation, even at 3 GPa. Calculation of G/Go using ab- 
initio density functional theory reveals that neither G nor the mobility, when properly averaged over 
the hole distribution, saturate at stresses lower than 3 GPa. The lack of saturation has important 
consequences for strained silicon technologies. 

PACS numbers: 73.50.Dn, 73.50.Gr, 73.63.Nm 



In semiconductors such as silicon, applied mechanical 
stress (X) modifies the band structure, principally chang- 
ing the carrier mobility and hence the conductance, G 
[l|-[3|- This phenomenon is known as piezoresistance. In 
n-type silicon, uniaxial compressive stress along the (100) 
crystal directions lifts the 6-fold degeneracy of the el- 
lipsoidal conduction band minima, resulting in a charge 
transfer from the four bands perpendicular to X into 
the two which are parallel with X. This charge transfer 
model [2] dictates an increase in the longitudinal con- 
ductivity effective mass (m^) that continues until all 
the electrons have transfered to the two parallel valleys 
at which point the resistance change saturates; it cor- 
rectly describes the piezoresistance of n-type silicon up 
to several hundred MPa [4[. The effect of stress on the 
valence band of silicon is more complex and cannot be de- 
scribed by a model of equivalent simplicity [5|. Instead, 
the degeneracy or near-degeneracy of the heavy (HH), 
light (LH) and split-off (SO) hole states requires the use 
of simplified analytical [3|, [6|, |7j or numerical techniques 
[8l4l2j that describe band warping and hole transfer ef- 
fects. For stresses up to several hundred MPa there is a 
consensus that the mobility (and hence G) change pri- 
marily due to a modification of ttiq 1: although some au- 
thors claim that changes to the inverse scattering rates 
(r m ) can be important [8|,LL3j. When using the band edge 
values of rrf GV these models are in reasonable agreement 
with piezoresistance data up to several hundred MPa. In 
the GPa stress range where there is as yet no experimen- 
tal data, a saturation of the piezoresistance is predicted 
once m* Gl ceases to change with X [ll|, LL3[, somewhat 
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reminiscent of the n-type (and germanium) cases 0, EH- 
In contrast, full band Monte Carlo methods [12| predict 
no such saturation. 

Coupled to these fundamental challenges and discrep- 
ancies in describing the piezoresistance of p-type silicon 
is the fact that GPa-range process induced stress has be- 
come technologically important since it can be used to re- 
duce tUq^ thereby improving device speeds and transcon- 
ductance gains of field-effect transistors [ll|, [l3|, [l5j . In- 
deed stressed silicon is one of the key elements of the 
ITRS roadmap [16|. The prediction of saturation in the 
hole mobility (i.e. in the piezoresistance) is thus an im- 
portant one since it will ultimately limit the increase in 
transistor performance afforded by strained-silicon tech- 
nology. 

Here we present piezoresistance data in n- and p-type 
silicon up to 3 GPa, well above previously experimen- 
tally accessible uniaxial stresses, and well beyond the 
predicted onset of saturation. While the piezoresistance 
saturation predicted by the charge transfer model holds 
in n-type silicon up to at least 2.4 GPa, in p-type silicon 
at even higher stresses no such saturation is observed, in 
stark contrast to theoretical predictions based on band 
edge effective masses. 

Experimentally, great care must be taken in order 
to reproducibly apply such large, mechanical stresses. 
To do so, a mechanical press was adapted to permit 
two-terminal conductance measurements perpendicularly 
across bulk semiconducting samples under uniaxial, com- 
pressive stress. The assembly depicted in Fig. [TJi ensures 
that the force is applied perpendicularly, by means of a 
flexible O-ring and a teflon sleeve that maintain paral- 
lelicity between the opposing steel chucks. The O-ring, 
being much softer than the sample, does not affect the ap- 
plied stress. Flat, mechanically polished tungsten plates 
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FIG. 1: (a) Cross-sectional view of the press, chuck assembly 
and sample, along the with the electrical connections for the 
conductance measurement. A compressive stress, X, is ap- 
plied parallel to the direction of current flow, (b) A series of 
photographs of a p-type (111) sample before each application 
of stress, showing the cleaving of the sides. 



are placed between the sample and the chucks to avoid in- 
dentation of the steel, as had occurred during initial test- 
ing of the assembly. Rigid, metal-covered non-conductors 
(ceramics and glass) were initially trialled in place of the 
tungsten but were found to reach their elastic limit and 
shatter at relatively low stress. The chucks themselves 
are electrically isolated from the press using glass plates 
to avoid short circuiting the sample with the press hous- 
ing. A load cell measures the applied force, which is 
recorded together with G. 

The samples were prepared from commercially avail- 
able double-side-polished n-type [(100), 5-10 D.cm, thick- 
ness 280 /im] and p-type [(110), > 1 l].cm, thickness 
400 /im] and [(111), 4-7 ft.cm, thickness 275 /im] sili- 
con wafers. Low doping densities were chosen to ensure 
that the measured two-terminal resistance is dominated 
by the sample itself rather than the leads. Ohmic con- 
tacts (p c = 10 -6 (].cm 2 ) are formed by a 100 nm deep 
thermally activated implant process on both faces of the 
wafer (doping density ~ 10 20 cm -3 ) and finished with 
annealed aluminum contacts. After processing, the sam- 
ples are cleaved with a diamond tipped scribe into pieces 
of size « 5 x 10 _7 m 2 (see Fig. [TJd), so that a pressure of 
1 GPa is obtained when about 500 N is applied. 

Before measurement, the surface area (S) of each sam- 
ple is estimated from an optical microscope image. In 
a first cycle, G is measured while the force is gradually 
applied to the sample until a discontinuous decrease in 
G is observed (point (a) in Fig. [2^). This jump is ac- 
companied by an audible crack and corresponds to the 
edges of the sample cleaving off in the press, thereby re- 
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FIG. 2: (a) Raw data obtained for X parallel to the (100) 
crystal direction in n-type silicon, showing discontinuous de- 
creases in conductance due to in-situ cleaving of the sample 
edges, (b) Data corrected for changes in sample surface area 
collapse onto a single curve that exhibits saturation above 1 
GPa. The fit of the charge transfer model to these curves 
(solid line) yields 'Eu — 12.1 eV and p±/p\\ = 2.9. The in- 
set shows the measured tti coefficient for this sample (black 
diamonds) and all other samples (open circles). 

ducing S and decreasing G (see Fig. [TJd). The reduction 
in S also leads to an instantaneous increase in X since 
the applied force remains unchanged. In the next cycle, 
the force is gradually ramped down to zero (without re- 
moving the sample) and then slowly ramped up until the 
sample cracks again. This process is repeated until the 
entire sample suddenly shatters into a fine powder, pre- 
sumably because the instantaneous increase in stress re- 
sulting from one cracking event causes a cascade of addi- 
tional cracking events. This typically occurs after about 
5 cycles. The resulting raw data is a series of curves re- 
lating G to the voltage output of the force transducer, 
as shown for a (lOO)-oriented n-type sample in Fig. [2^. 
These curves cannot be directly compared because S is 
different for each one, so that the relationship between 
X and the transducer output is not fixed. Furthermore, 
the zero-stress conductance (Go) is reduced as S is de- 
creased. S is calculated for each cycle by multiplying the 
initial value by the ratio of the present and initial val- 
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lies of Go- In this way, the measured force and S can 
be used to determine X for each cycle. The conductance 
during each sweep is also normalized to Go, allowing a 
direct comparison of the piezoresistance for different sur- 
face areas. The corrected G/Go data are plotted against 
X in Fig. Eb. The curves lie very close to one another 
and saturation occurs at « 1 GPa. 
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FIG. 3: Corrected data obtained along the (a) (110) and (b) 
(111) directions in p-type silicon (black symbols). No satura- 
tion is observed in either case. The insets show the measured 
tti coefficients for all samples. In both cases the result of the 
ab-initio calculation is shown when concentration weighted 
band edge m* Gl for the HH and LH bands are used (lines with 
open circles), and when Eq. [T]is used to average over the LH 
and HH hole distributions (lines with closed circles). 

The data in Fig. [2]o is fitted using Aubrey's adapta- 
tion [4] of the charge transfer model [2] which has two 
free parameters; S u , the shear deformation potential and 
, the mobility ratio of the ellipsoidal valleys per- 
pendicular and parallel to X. The solid curve in Fig. 
[2]o corresponds to this fit and, over all the n-type sam- 
ples tested, E u = 12.6 ± 3.5 eV and fJL±/fJL\\ = 2.6 ± 0.4 
is found. These values correspond well with those cal- 
culated and measured elsewhere [H, [3, 0] , indicating that 
the experimental method and the conductance correction 
technique described above are reasonable. Furthermore, 
the slopes of the G/Go curves at low stress yield the lon- 
gitudinal piezoresistance coefficient, 717, from which it is 
found that m = dG/dX x 1/G = -82 ± 15 x lO" 11 



Pa -1 which is slightly lower than the generally accepted 
value [1]. While the difference is not excessive, espe- 
cially when compared to the know experimental scatter 
in reported values for the [100] crystal direction pj), it 
may be the result of an unintentionally applied trans- 
verse stress. Given the small aspect ratio of the sample 
it is conceivable that transverse stress cannot relax over 
the sample thickness. Since the transverse piezoresis- 
tance coefficients have the opposite sign to 717, this may 
systematically result in slightly lower apparent 717. 

Figure [3^ and b show the data for p-type silicon with 
X || (110) and X || (111) crystal directions respec- 
tively. For I < 1 GPa the slopes of these curves yield 
717 = 140 ± 21 x 10" 11 Pa -1 for the (110) direction and 
717 = 159 ± 37 x 10" 11 Pa -1 for the (111) direction (see 
insets in Fig. [3^i and Fig. [3}d) which is a factor of two 
higher than the generally accepted low-stress values [l|. 
In attempting to understand this difference it is impor- 
tant to note that in Smith's original work X is a ten- 
sile stress between 1 and 10 MPa i.e. smaller than the 
stress increment in the GPa-adapted apparatus described 
here. Moreover there is some experimental [l8|, [l9| and 
theoretical [l(| G3 evidence in the literature that the 
p-type piezoresistance along the [110] direction is larger 
for compressive stresses than for tensile stresses. One 
experimental observation reports values of 717 under 600 
MPa of compressive stress that are 1.5 times larger than 
Smith's value [l8j . This is in excellent agreement with 
a theory [12] which predicts increases of a factor of 2 in 
717 at higher stresses i.e. similar to the values measured 
here. The intermediate stress (X < 600 MPa) values of 
717 reported here (at least for the [110] crystal direction) 
should probably not therefore be compared directly with 
Smith's values, but rather with those of Shifren [18| and 
Tsang [19|. The remaining difference may be the result 
of a lack of parallelicity in the tungsten contacts with 
the sample surfaces. In this case X will only be applied 
locally to parts of the sample, and will be larger than its 
averaged value. This will therefore tend to systematically 
underestimate X and overestimate 717 . A final (and most 
important) aspect of the data is the lack of saturation at 
high stress for both crystal orientations (see dotted hor- 
izontal line in Fig. [3^ for the predicted saturation level, 
[U, 0). At X « 3 GPa, G/Go ~ 6.5, well above the 
predicted limit of about 4.5. If X is systematically un- 
derestimated due to contact parallelicity problems then 
this result is even more surprising since X is at least 3 
GPa. 

In order to investigate the discrepancy between the 
predicted saturation in G/Go and the experimental ob- 
servations, ab-initio density functional calculations of the 
full Brillouin zone band structure were performed up to 
3 GPa (oriented along the (100), (110) and (111) crystal 
directions) using ABINIT [20| with the PAW method to 
account for spin-orbit coupling (2l|. The calculation is 
performed by imposing a target stress tensor and then 
performing a structural optimization which yields the 
deformed lattice parameters. These are then used to 
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calculate the band structure and the density of states. 
The tetrahedron integration method is modified in order 
to obtain the partial density of states for the i th band, 
Di(E). In the conduction band with X ||< 100 >, G/Go 
is in excellent agreement with the charge transfer model 
shown in Fig. EJd (ab-initio result not shown). 

In the following discussion of the valence band with 
X ||< 110 > or X || < 111 >, the total hole concentration 
is assumed to be a stress-independent p = 10 16 cm -3 
which is imposed using the calculated Di{E) and fixing 
the Fermi energy as necessary. The GW correction [22| is 
used to confirm that the intrinsic concentration change 
due to the applied stress is indeed negligible compared 
to p. It is also noted that since the semiconductor is in 
the non-degenerate limit, the calculated values of G/Go 
depend only very weakly on the exact value of p. 

For X || < 110 > the band edge calculation of 
for the HH and LH bands, weighted by the respective 
concentrations, results in a saturation of G/Go at « 3.5 
(see lines with open circle symbols in Fig. [3^i) in qual- 
itative agreement with the behavior predicted elsewhere 
[l3| (see dotted line in Fig. [3^). This clearly does not 
describe the experimental data. In order to better do 
so, the calculated dispersion relations for the two bands 
which contribute significantly to the transport (HH and 
LH) are used to find the velocity as a function of energy, 
Vi(E) oc dE/dk\i for the i th band. The partial conduc- 
tance of the i th band is then proportional to Koi, the 
usual zeroth order moment of the solution to the Boltz- 
mann transport equation [23|: 

K 0i oc J $(E)T mi (E)^Di(E)dE. (1) 

Here df°/dE is the Fermi window function with f° the 
equilibrium Fermi-Dirac distribution function. In order 
to properly evaluate Eq. [TJ the energy dependence of 
the relaxation time, r m i(E), must be calculated. Here 
an oversimplified but useful approach is taken: it is as- 
sumed that momentum relaxation is determined by in- 
traband optical phonon scattering. Given the angular 
dependence of the interband and intraband overlap fac- 
tors for the scattering rates [24J , the error induced by this 
approximation is probably not too large. In this case [23| 
T m i(E) oc l/Di(E) and Ko% no longer depends on Di(E). 
G is then calculated according to G oc J^-Kio- The re- 
sulting dependence of G/Go on X (shown in Fig. [3^, lines 
with solid circles) exhibits no saturation and is surpris- 
ingly close to the measured data. It is noted that a full 
band Monte Carlo simulation [12] yields similar results 
that are in excellent agreement with the experimental 
data presented here. 

Fig. 2^ details the calculated changes in band struc- 
ture oriented along a (110) direction, parallel to X. In 
this case the LH (HH) band is represented by open 
(closed) circles for three values of X: 0, 1.5 and 3 GPa. 
Stress increases the energy of the HH band maximum at 
k = and results in a LH/HH anti-crossing for k ^ 0. 
For X < 1 GPa, the anti-crossing is close to the zone 
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FIG. 4: (a) The LH (open circles) and HH (closed cir- 
cles) band structure calculated along (110) parallel to X for 
three stress levels indicated by the circle size: (small), 1.5 
(medium) and 3 GPa (large). LH/HH anti-crossing plays an 
important role in determining ensemble averaged G for X < 1 
GPa where the integrand of Eq. [T]is non- negligible (b). At 
high stresses, holes progressively occupy only the HH result- 
ing in a gradual roll-off in G without hard saturation. 



center and significantly modifies the ensemble averaged 
effective transport mass (and thus v(E)) for both bands. 
This is important for the overall value of G since the in- 
tegrand of Eq. [1] (see Fig. HJd) is non-negligible close to 
the anti-crossing. At higher stresses the HH/LH split- 
ting increases and the anti-crossing moves out to larger k 
where the integrand becomes smaller due to the finite size 
of the Fermi window function. This reflects the physical 
fact that as X increases, holes are progressively restricted 
to the HH band at energies well above the anti-crossing. 
It is also the reason for the lack of hard saturation in 
G/Go since the thermal distribution of holes does not 
abruptly fall to zero close to the band edge, but is rather 
smeared out over the band dispersion. In this very sim- 
ple description a saturation in G/Go is still expected as 
can be seen in the slight roll-off in the calculated curve in 
Fig. [3^, although hard saturation will presumably only 
occur for X ^> 3 GPa (note also that indeed a weak onset 
of saturation is sometimes observed as in Fig. [3^ (open 
diamonds)). 

For X || (111) the behaviour is qualitatively similar to 
the (110) direction. A saturation in G/Go is predicted 
(see lines with open circles in Fig. (3}d) when accounting 
only for the concentration weighted (HH and LH) band 
edge values of m* Gl whereas the ensemble average shows 
no saturation whatsoever. Again, the application of Eq. 
[1] yields a surprisingly good agreement with the experi- 
mental data (see lines with solid circles in Fig. [3]o) . 

The concrete experimental evidence presented here 
along with a simple theoretical interpretation strongly 
suggest that up to 3 GPa at least, the variation of G/Go 
with X is due entirely to the ensemble averaged mobil- 
ity variation, and not to a change in p. Moreover, it 
also suggests the mobility variation arises from both a 
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change in m* Gl due to charge transfer and band warp- 
ing, and to an increase in r m with compressive stress. 
The qualitative agreement between experiment and the- 
ory tentatively indicates that the interband scattering 
contribution to r m is negligible. These conclusions have 
potentially important consequences for p-type strained- 
silicon devices, particularly in the technologically impor- 
tant (110) directions. They indicate that much larger 
hole mobilities than previously envisaged are possible at 
process-induced stresses higher that those which are cur- 
rently accessible in p-type transistor channels, suggesting 
that further investment in this technology may yield div- 



idends. 
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